THE GEOMETRIC STABILITY OF VORONOI DIAGRAMS 
WITH RESPECT TO SMALL CHANGES OF THE SITES 



DANIEL REEM 



Abstract. Voronoi diagrams appear in many areas in science and technol- 
ogy and have numerous apphcations. They have been the subject of extensive 
investigation during the last decades. Roughly speaking, they are a certain 
decomposition of a given space into cells, induced by a distance function and 
by a tuple of subsets called the generators or the sites. Consider the follow- 
ing question: does a small change of the sites, e.g., of their position or shape, 
yield a small change in the corresponding Voronoi cells? This question is by 
all means natural and fundamental, since in practice one approximates the sites 
either because of inexact information about them, because of inevitable numer- 
ical errors in their representation, for simplification purposes and so on, and it 
is important to know whether the resulting Voronoi cells approximate the real 
ones well. The traditional approach to Voronoi diagrams, and, in particular, to 
(variants of) this question, is combinatorial. However, it seems that there has 
been a very limited discussion in the geometric sense (the shape of the cells), 
mainly an intuitive one, without proofs, in Euclidean spaces. We formalize this 
question precisely, and then show that the answer is positive in the case of M.'^, 
or, more generally, in (possibly infinite dimensional) uniformly convex normed 
spaces, assuming there is a common positive lower bound on the distance be- 
tween the sites. Explicit bounds arc given, and we allow infinitely many sites of 
a general form. The relevance of this result is illustrated using several pictures 
and many real-world and theoretical examples and counterexamples. 
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1. Introduction 

1.1. Background. The Voronoi diagram (the Voronoi tessellation, the Voronoi 
decomposition, the Dirichlet tessellation) is one of the basic structures in compu- 
tational geometry. Roughly speaking, it is a certain decomposition of a given space 
X into cells, induced by a distance function and by a tuple of subsets {Pk)k£K, 
called the generators or the sites. More precisely, the Voronoi cell Rk associated 
with the site Pk is the set of all the points in X whose distance to Pk is not greater 
than their distance to the union of the other sites Pj. 

Voronoi diagrams appear in a huge number of fields in science and technology 
and have many applications. They have been the subject of research for at least 
160 years, starting formally with L. Dirichlet [17] and G. Voronoi [58], and of 
extensive research during the last 40 years. For several well written surveys on 
Voronoi diagrams which contain extensive bibliographies and many applications, 
see [7], [8], [44], and [23]. 

Consider the following question: 

Question 1.1. Does a small change of the sites, e.g., of their position or shape, 
yield a small change in the corresponding Voronoi cells? 

This question is by all means natural, because in practice, no matter which al- 
gorithm is being used for the computation of the Voronoi cells, one approximates 
the sites either because of lack of exact information about them, because of in- 
evitable numerical errors occurring when a site is represented in an analog or a 
digital manner, for simplification purposes and so on, and it is important to know 
whether the resulting Voronoi cells approximate well the real ones. 

For instance, consider the Voronoi diagram whose sites are either shops (or large 
shopping centers), antennas, or other facilities in some city/district such as post 
offices. See Figures 1-2. 

Each Voronoi cell is the domain of influence of its site and it can be used for 
various purposes, among them estimating the number of potential costumers [57] 
or understanding the spreading patterns of mobile phone viruses [60]. In reality, 
each site has a somewhat vague shape, and its real location is not known exactly. 
However, to simplify matters we regard each site as a point (or a finite collection 
of points if we consider firms of shops) located more or less near the real location. 
As a result, the resulting cells only approximate the real ones, but we hope that 
the approximation will be good in the geometric sense, i.e., that the shapes of the 
corresponding real and approximate cells will be almost the same. (See Section 
4 for many additional examples, including ones with infinitely many sites or in 
higher /infinite dimensional spaces.) As the counterexamples in Section 6 show, it 
is definitely not obvious that this is the case. 

A similar question to Question 1.1 can be asked regarding any geometric struc- 
ture/algorithm, and, in our opinion, it is a fundamental question which is analo- 
gous to the question about the stability of the solution of a differential equation 
with respect to small changes in the initial conditions. 
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Figure 1. 10 shopping centers 
(or post offices) in a flat city. 
Each shopping center is repre- 
sented by a point. 




Figure 2. In reahty each shop- 
ping center /post office is not a 
point and its location is approx- 
imated. The combinatorial struc- 
ture is somewhat different and the 
Voronoi cells are not exactly poly- 
gons, but still, their shapes are al- 
most the same as in Figure 1. 



The traditional approach to Voronoi diagrams, and, in particular, to (variants 
of) Question 1.1, is combinatorial. For instance, as already mentioned in Auren- 
hammer [7, p. 366], the combinatorial structure of Voronoi diagrams (in the case 
of the Euclidean distance with point sites), i.e., the structure of vertices, edges and 
so on, is not stable under continuous motion of the sites, but it is stable "most of 
the time" . A more extensive discussion about this issue, still with point sites but 
possibly in higher dimensions, can be found in Weller [61], Vyalyi et al. [59], and 
Albers et al. [2]. 

However, it seems that this question, in the geometric sense, has been raised 
or discussed only rarely in the context of Voronoi diagrams. In fact, after a com- 
prehensive survey of the literature about Voronoi diagrams, we have found only 
very few places that have a very brief, particular, and intuitive discussion which 
is somewhat related to this question. The following cases were mentioned: the 
Euclidean plane with finitely many sites [32], the Euclidean plane with finitely 
many point sites [7, p. 366], and the ci-dimensional Euclidean space with finitely 
many point sites [2] (see also some of the references therein). It was claimed there 
without proofs and exact definitions that the Voronoi cells have a continuity prop- 
erty: a small change in the position or the shape of the sites yields a small change 
in the corresponding Voronoi cells. 

Another continuity property was discussed by Groemer [25] in the context of 
the geometry of numbers. He considered Voronoi diagrams generated by a lattice 
of points in a dimensional Euclidean space, and proved that if a sequence of 
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lattices converges to a certain lattice (meaning that the basis elements which 
generate the lattices converge with respect to the Euclidean distance to the basis 
which generates the limit lattice), then the corresponding Voronoi cells of the 
origin converge, with respect to the Hausdorff distance, to the cell of the origin of 
the limit lattice. His result is, in a sense and in a very particular stabihty 
result, but it definitely does not answer Question 1.1 (which, actually, was not 
asked at all in [25]) for several reasons: first, usually the sites or the perturbed 
ones do not form a (infinite) lattice. Second, in many cases they are not points 
(singletons). Third, a site is usually different from the perturbed site (in [25] the 
discussed sites equal {0}). In this connection, we also note that Groemer's proof 
is very restricted to the above setting and it uses arguments based on compactness 
and no explicit bounds are given. 

It is quite common in the computational geometry literature to assume "ideal 
conditions", say infinite precision in the computation, exact and simple input, and 
so on. These conditions are somewhat non-realistic. Issues related to the stability 
of geometric structures under small perturbations of their building blocks (not nec- 
essarily the geometric stability) are not so common in the literature, but they can 
be found in several places, e.g., in [1, 3, 6, 9, 11, 12, 15, 21, 30, 33, 40, 41, 50, 55]. 
However, in many of the above places the discussion has combinatorial character- 
istics and there are several restrictive assumptions: for instance, the underlying 
setting is usually a finite dimensional space (in many cases only or R'^), with 
the Euclidean distance, and with finitely many objects of a specific form (merely 
points in many cases). In addition, the methods are restricted to this setting. In 
contrast, the infinite dimensional case or the case of (possibly infinitely many) 
general objects or general norms have never been considered. 

1.2. Contribution of this paper. We discuss the question of stability of Voronoi 
diagrams with respect to small changes of the corresponding sites. We first formal- 
ize this question precisely, and then show that the answer is positive in the case 
of M'^, or, more generally, in the case of (possibly infinite dimensional) uniformly 
convex normed spaces, assuming there is a common positive lower bound on the 
distance between the sites. Explicit bounds are presented, and we allow infinitely 
many sites of a general form. We also present several counterexamples which show 
that the assumptions formulated in the main result are crucial. We illustrate the 
relevance of this result using several pictures and many real-world and theoretical 
examples and counterexamples. To the best of our knowledge, the main result 
and the approach used for deriving it are new. Two of our main tools are: a 
new representation theorem which characterizes the Voronoi cells as a collection 
of line segments and a new geometric lemma which provides an explicit geometric 
estimate. 

1.3. The structure of the paper. In Section 2 we present the basic definitions 
and notations. Exact formulation of Question 1.1 and informal description of the 
main result are given in Section 3. The relevance of the main result is illustrated 
using many theoretical and real-world examples in Section 4. The main result is 
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presented in Section 5, and we discuss briefly some aspects related to its proof. 
In Section 6 we present several interesting counterexamples showing that the as- 
sumptions imposed in the main result are crucial. We end the paper in Section 7 
with several concluding remarks. Since the proof of the main result is quite long 
and technical, and because the main goal of this paper is to introduce the issue 
and to discuss it in a qualitative manner, rather than going deep into technical 
details, proofs were omitted from the main body of the text. Full proofs can be 
found in the appendix (Section 8) and a preliminary version in [47]. 

2. Notation and basic definitions 

In this section we present our notation and basic definitions. In the main dis- 
cussion we consider a closed and convex set X 7^ in some uniformly convex 
normed space {X, | ■ |) (see Definition 2.3 below), real or complex, finite or infi- 
nite dimensional. The induced metric is d{x,y) = \x — y\. We assume that X 
is not a singleton, for otherwise everything is trivial. We denote by [p, x] and 
\p,x) the closed and half open line segments connecting p and x, i.e., the sets 
{p + t{x — p) : t E [0,1]} and {p + t{x —p) : t E [0, 1)} respectively. The (possibly 
empty) boundary of X with respect to the affine hull spanned by X is denoted by 
dX . The open ball with center x E X and radius r > is denoted by B{x, r). 

Definition 2.1. Given two nonempty subsets P, A (1 X , the dominance region 
dom(P, A) of P with respect to A is the set of all x E X whose distance to P is 
not greater than their distance to A, i.e., 

dom(P, A) = {xEX : d{x, P) < d{x, A)}. 

Here d{x,A) = ini{d{x,a) : a E A} and in general we denote ^(741,^2) = 
inf{(i(ai, 02) : Oi E Ai, a2 E A2} for any nonempty subsets Ai,A2. 

Definition 2.2. Let K be a set of at least 2 elements (indices), possibly infinite. 
Given a tuple {Pk)k&K of nonempty subsets Pk C X, called the generators or the 
sites, the Voronoi diagram induced by this tuple is the tuple {Rk)k£K of non-empty 
subsets Rk C X, such that for all k E K , 

Rk = dom(Pfc, IJP,) = {xeX: rf(x, Pk) < d{x, P,) Vj E K, j ^ k}. 

In other words, the Voronoi cell Rk associated with the site Pk is the set of all 
X E X whose distance to Pk is not greater than their distance to the union of the 
other sites Pj . 

In general, the Voronoi diagram induces a decomposition of X into its Voronoi 
cells and the rest. If K is finite, then the union of the cells is the whole space. 
However, if K is infinite, then there may be a "neutral cell" : for example, if X is 
the Euclidean plane, K = N = {1,2,3, . . .} and P^ = M x {1/k}, then no point 
in the lower half-plane M x (—00, 0] belongs to any Voronoi cell. In the above 
definition and the rest of the paper we ignore the neutral cell. 

We now recall the definition of strictly and uniformly convex spaces. 
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Definition 2.3. A normed space {X, \ ■ \) is said to be strictly convex if for all 
x,y E X satisfying \x\ = \y\ = 1 and x ^ y, the inequality \{x + ?/)/2| < 1 holds. 
[X, \ ■ \) is said to he uniformly convex if for any e G (0,2] there exists 5 G (0, 1] 
such that for all x,y E X , if \x\ = \y\ = 1 and \x — y\ > e, then \ {x + y)/2\ < 1 — 5. 

Roughly speaking, if tlie space is uniformly convex, then for any e > there 
exists a uniform positive lower bound on how deep the midpoint between any two 
unit vectors must penetrate the unit ball, assuming the distance between them is 
at least e. In general normed spaces the penetration is not necessarily positive, 
since the unit sphere may contain line segments. M? with the max norm | ■ |oo is a 
typical example for this. A uniformly convex space is always strictly convex, and 
if it is also finite dimensional, then the converse is true too. The m-dimensional 
Euclidean space M*", or more generally, inner product spaces, the sequence spaces 
£p, the Lebesgue spaces LpiVL) {1 < p < oo), and a uniformly convex product of 
a finite number of uniformly convex spaces, are all examples of uniformly convex 
spaces. See Clarkson [13] and, for instance, Goebel-Reich [22] and Lindenstrauss- 
Tzafriri [38] for more information about uniformly convex spaces. 

From the definition of uniformly convex spaces we can obtain a function which 
assigns to the given e a corresponding value 5(e). There are several ways to obtain 
such a function, but for our purpose we only need 5 to be increasing, and to satisfy 
(5(0) = and 5(e) > for any e G (0,2]. One choice, which is not necessarily the 
most convenient one, is the modulus of convexity, which is the function 5 : [0, 2] — > 
[0, 1] defined by 

(5(e) = inf{l — \{x + y)/2| : \x — y\> e, \x\ = \y\ = 1}. 

For specific spaces we can take more convenient functions. For instance, for the 
spaces Lp{Q) or £p , 1 < p < oo, we can take 

5(e) = 1 - (1 - (e/2)P)^/^ forp>2, 

5(e) = 1 - (1 - (e/2)'')^/'' , for 1 < p < 2 and i + i = 1. 

We finish this section with the definition of the Hausdorff distance, a definition 
which is essential for the rest of the paper. 

Definition 2.4. Let {X,d) be a metric space. Given two nonempty sets Ai,A2 C 
X , the Hausdorff distance between them is defined by 

D{Ai,A2) = max{ sup d{ai,A2), sup d{a2,Ai)}. 

aiGAi a2&A2 

Note that the Hausdorff distance D{Ai., A2) is definitely different from the usual 
distance d{Ai,A2) = inf{(i(ai, 02) : ai G Ai, 02 G A2}. As a matter of fact, 
D{Ai,A2) < e if and only if d{ai,A2) < e for any ai G Ai, and d{a2,Ai) < e for 
any 02 G A2. In addition, if D{Ai,A2) < e, then for any ai G Ai there exists 
02 G A2 such that d{ai, 02) < e, and for any 62 ^ ^2 there exists bi G Ai such that 
d{b2,bi) < e. These properties explain why the Hausdorff distance is the natural 
distance to be used when discussing approximation and stability in the context of 
sets: suppose that our resolution is at most r, i.e., we are not able to distinguish 
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between two points whose distance is at most some given positive number r. If it 
is known that D{Ai, A2) < r, then we cannot distinguish between the sets Ai and 
A2, at least not by inspections based only on distance measurements. As a result 
of the above discussion, the intuitive phrase "two sets have almost the same shape" 
can be formulated precisely: the Hausdorff distance between the sets is smaller 
than some given positive parameter (note that a set and a rigid transformation of 
it usually have different shapes). 

3. Exact formulation of the main question and informal 
formulation of the main result 

The exact formulation of Question 1.1 is based on the concept of Hausdorff 
distance for reasons which were explained at the end of the previous section. 

Question 3.1. Suppose that {Pk)k£K is a tuple of non-empty sets in X. Let 
{Rk)k£K be the corresponding Voronoi diagram. Is it true that a small change 
of the sites yields a small change in the corresponding Voronoi cells, where both 
changes are measured with respect to the Hausdorff distance? More precisely, is 
it true that for any e > there exists A > such that for any tuple {Pl,)k£K, 
the condition D{Pk, PI.) < A for each k E K implies that D{Rk, R'l^) < e for each 
k E K, where {R'f^jkeK is the Voronoi diagram of {Pl,)keK'^ 

The main result (Theorem 5.1) says that the answer is positive. Here is an 
informal description of it: 

Answer 3.2. Suppose that the underlying subset X is a closed and convex set 
of a (possibly infinite dimensional) uniformly convex normed space X. Suppose 
that a certain boundedness condition on the distance between points in X and the 
sites holds, e.g., when X is bounded or when the sites form a (distorted) lattice. 
If there is a common positive lower bound on the distance between the sites, and 
the distance to each of them is attained, then indeed a small enough change of the 
(possibly infinitely many) sites yields a small change of the corresponding Voronoi 
cells, where both changes are measured with respect to the Hausdorff distance; in 
other words, the shapes of the real cells and the corresponding perturbed ones are 
almost the same. Moreover, explicit bounds on the changes can be derived and 
they hold simultaneously for all the cells. There are counterexamples which show 
that the assumptions imposed above are crucial. 

The condition that the distance to a site is attained holds, e.g., when the site 
is either a closed subset contained in a (translation of a) finite dimensional space, 
or a compact set, or a convex and closed subset in a uniformly convex Banach 
space. The sites can always be assumed to be closed, since the distance and the 
Hausdorff distance preserve their values when the involved subsets are replaced 
by their closures. The "certain boundedness condition on the distance between 
points in X and the sites" is a somewhat technical condition expressed in (2) (see 
also Remark 5.2). 
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4. The relevance of the main result 

In Section 1 we explained why Question 1.1 is natural and fundamental, and 
mentioned the real- world example of a Voronoi diagram induced by shops/cellular 
antennas. The goal of this section is to illustrate further the relevance of the 
main result using a (far from being exhaustive) list of real-world and theoretical 
exampls. In these examples the shape or the position of the real sites are obviously 
approximated, and the main result (Theorem 5.1) ensures that the approximate 
Voronoi cells and the real ones have almost the same shape, and no unpleasant 
phenomenon such as the one described in Figures 5-8 can occur. 

One example is in molecular biology for modeling the proteins structure (Richards 
[48], Kim et al. [34], Bourquard et al. [10]), where the sites are either the 
atoms of a protein or special selected points in the amino acids and they are 
approximated by spheres/points. Another example is related to collision detec- 
tion and robot motion (Goralski-Gold [24], Schwartz et al. [51]), where the sites 
are the (static or dynamic) obstacles located in an environment in which a vehi- 
cle/airplane/ship/robot/satellite should move. A third example is in solid state 
physics (Ashcroft-Mermin [5]; here the common terms are "the first Brillouin zone" 
or "the Wigner-Seitz cell" instead of "the Voronoi cell"), where the sites are in- 
finitely many point atoms in a (roughly) periodic structure which represents a 
crystal. A fourth example is in material engineering (Li-Ghosh [36]), where the 
sites are cracks in a material. 

A fifth example is in numerical simulations of various dynamical phenomena, 
e.g., gas, fluid or urban dynamics (Slotterback et al. [53], Mostafavi et al. [43]). 
Here the sites are certain points/shapes taken from the sampled data of the simu- 
lated phenomena, and the cells help to cluster and analyze the data continuously. 
A sixth example is in astrophysics (Springel et al. [54]) where the (point) sites are 
actually huge regions in the universe (of diameter equals to hundreds of light years) 
used in simulations performed for understanding the behavior of (dark) matter in 
the universe. A seventh example is in image processing and computer graphics, 
where the sites are either certain important features/parts in an image (Tagare et 
al. [56], Dobashi et al. [18], Sabha-Dutre [49]) used for processing/analyzing it, 
or they are points form a useful configuration such as (an approximate) centroidal 
Voronoi diagram (CVD) which induces cells having good shapes (Du et al. [19], 
Liu et al. [39], Faustino-Figueiredo [20]). 

An eighth example is in computational geometry, and it is actually a large 
collection of familiar problems in this field where Voronoi cells appear and being 
used, possibly indirectly: (approximate) nearest neighbor searching/the post office 
problem, cluster analysis, (approximate) closest pairs, motion planning, finding 
(approximate) minimum spanning trees, finding good triangulations, and so on. 
See, e.g., Aurenhammer [7], Aurenhammer-Klein [8], Clarkson [14], Indyk [31], 
and Okabe et al. [44]. Here the sites are either points or other shapes, and the 
space is usually R"" with some norm. In some of the above problems our stability 
result is clearly related because of the analysis being used (e.g., cluster analysis) or 
because the position/shapes of the sites are approximated (e.g., motion planning. 
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the post office problem). However, it may be related also because in many of 
the previous problems the difficulty/complexity of finding an exact solution is too 
high, so one is forced to use approximate algorithms, to impose a general position 
assumption, and so on. Now, by perturbing slightly the sites (moving points, 
replacing a non-smooth curve by an approximating smooth one, etc.,) one may 
obtain much simpler and tractable configurations and now, using the geometric 
stability of the Voronoi cells, one may estimate how well the obtained solution 
approximates the best/real one. 

As for a theoretical example of a different nature, we mention Kopecka et al. 
[35] in which the stability results described here have been used, in particular, 
for proving the existence of a zone diagram (a concept which was first introduced 
by Asano et al. [4] in the Euclidean plane with point sites) of finitely many 
compact sites which are strictly contained in a (large) compact and convex subset 
of a uniformly convex space, and also for proving interesting properties of Voronoi 
cells there. 

Another example is for the infinite dimensional Hilbert space L2{I) for some 
/ (perhaps an interval or a finite dimensional space): functions in it are used in 
signal processing and in many other areas in science and technology. In practice 
the signals (functions) can be distorted, e.g., because of noise, and in addition, 
they are approximated by finite series (e.g., finite Fourier series) or integrals (e.g., 
Fourier transform). Given several signals, the (approximate) Voronoi cell of a 
given signal may help, at least in theory, to cluster or analyze data related to the 
sites. Such an analysis can be done also when the signal is considered as a point in 
a finite dimensional space of a high dimension. See, for instance, Conway-Sloane 
[16, pp. 66-69, 451-477] (coding) and Shannon [52] (communication) for a closely 
related discussion (in [52] Voronoi diagrams are definitely used in various places, 
but without their explicit name). 

We mention several additional examples related to our stability result, some- 
times in a somewhat unexpected way. For instance, Voronoi diagrams of infinitely 
many sites generated by a Poisson process (Okabe et al. [44, pp. 39, 291-410]), 
Voronoi diagrams of atom nuclei used for the mathematical analysis of stability 
phenomena in matter (Lieb-Yau [37]), Voronoi diagrams of infinitely many lattice 
points in a multi-dimensional Euclidean space which appear in the original works 
of Dirichlet [17] and Voronoi [58] (see also Groemer [25] and Gruber-Lekkerkerker 
[26] regarding the geometry of numbers and quadratic forms; Groemer used his 
stability result for deriving the Mahler compactness theorem [42]), and packing 
problems such as the Kepler conjecture and the Dodecahedral conjecture (Hales 
[27], [28], Hales-McLaughlin [29]; because of continuity arguments needed in the 
proof) or those described in Conway-Sloane [16]. 

5. The main result and some aspects related to its proof 

In this section we formulate the main result and discuss briefly issues related to 
its proof. See also the remarks after Theorem 5.1 for several relevant clarifications. 
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Figure 3. Illustration of Theo- 
rem 5.1: five sites in a square 
in (M^, ip) where the parameter is 
p = 3.14159. 




Figure 4. The sites have been 
slightly perturbed: the two points 
have merged, the "sin" has 
shrunk, and so on. The cells have 
been slightly perturbed. 



Theorem 5.1. Let (X, | ■ |) he a uniformly convex normed space. Let X ^ X be 
closed and convex. Let {Pk)keK, {P'k)k£K be two given tuples of nonempty subsets 
of X with the property that the distance between each x E X and each Pk, PI. is 
attained. For each k E K let = Ujyfc-^i' ^fc ~ Uf 
following conditions hold: 



Suppose that the 



(1) 



r/ := mi{d{Pk, P,) ■.j,ke K,j ^ fc} > 0, 



(2) 3p G (0, oo) such that for all k E K and for all x G X 

the open ball B{x,p) intersects A^. 

For each k E K let = dom(Pfc, A^), i?^ = dom(P^,yl'^) be, respectively, the 
Voronoi cells associated with the original site Pk and the perturbed one P^. Then 
for each e G (0, 77/6) there exists A > such that if D{Pk, P^) < A for each k E K, 
then D{Rk, R'^) < e for each k E K . 

See Figures 3, 4 for an illustration. The pictures were produced using the 
algorithm described in [46]. 

Remark 5.2. The assumption mentioned in (2) may seem somewhat complicated 
at a first glance, but it actually expresses a certain uniform boundedness condition 
on the distance between any point in X to its neighbor sites. No matter which 
point X E X and which site Pk are chosen, the distance between x and the collection 
of other sites Pj, j ^ k cannot be arbitrary large. A sufficient condition for it to 
hold is when a uniform bound on the diameter of the cells (including the neutral 
one, if it is nonempty) is known in advance, e.g., when X is bounded or when the 
sites form a (distorted) lattice. But (2) can hold even if the cells are not bounded, 
e.g., when the setting is the Euclidean plane and P^ = M x {k} where k runs over 
all integers. 
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Remark 5.3. In general, we have A = O(e^). However, if there is a positive lower 
bound on the distance between the sites and the boundary of X (relative to the 
affine hull spanned by X), i.e., if the sites are strictly contained in the (relative) 
interior of X, then actually the better estimate A = 0(e) holds. The constants 
inside the big O can be described explicitly: when A = Ce^ we can take 



In the second case, in addition to e < 77/6, the inequality e < 8 ■ d(\j^^j^ Pk.dX) 
should be satisfied too. 

The proof of Theorem 5.1 is quite long and technical, and hence it is given in the 
appendix. Despite this, we want to say a few words about the proof and some of 
the difficulties which arise along the way. First, as the counterexamples mentioned 
in Section 6 show, one must take into account all the assumptions mentioned in 
the formulation of the theorem. 

Second, in order to arrive to the generality described in the theorem, one is forced 
to avoid many familiar arguments used in computational geometry and elsewhere, 
such as Euclidean arguments (standard angles, trigonometric functions, normals, 
etc.), arguments based on lower envelopes and algebraic arguments (since the in- 
tersections between the graphs generating the lower envelope may be complicated 
and since the boundaries of the cells may not be algebraic), arguments based 
on continuous motion of points, arguments based on finite dimensional properties 
such as compactness (since in infinite dimensional spaces closed and bounded balls 
are not compact), arguments based on finiteness (since we allow infinitely many 
sites and sites consist of infinitely many points) and so on. Our way to overcome 
these difficulties is to use a new representation theorem for dominance regions as 
a collection of line segments (Theorem 5.4 below) and a new geometric lemma 
(Lemma 5.5 below) which enables us to arrive to the explicit bounds mentioned 
in the theorem. As a matter of fact, we are not aware of any other way to obtain 
these explicit bounds even in a square in the Euclidean plane with point sites. 
These tools also allow us to overcome the difficulty of a potential infinite accumu- 
lated error due to the possibility of infinitely many sites/sites with infinitely many 
points/ infinite dimension. 

Theorem 5.4. Let X he a closed and convex subset of a normed space (X, | ■ |), 
and let P, A C X he nonempty. Suppose that for all x & X the distance hetween 
X and P is attained. Then dom(P, A) is a union of line segments starting at the 
points of P. More precisely, given p & P and a unit vector 9, let 



T{9,p) = sup{t G [0,00) : p + t9 e X and d{p + t9,p) < d{p + t9,A)}. 




and when A = Ce we can take 
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Then 



dom{P,A)= [j [j[p,p + T{9,p)9]. 



peP \e\=i 



WhenT{6,p) = 00, the notation [p,p+T {6, p)6] means the ray {p+t9 : t E [0, 00)}. 

Lemma 5.5. Let (X, | ■ |) he a uniformly convex normed space and let A (1 X be 
nonempty. Suppose that y,p E X satisfy d{y,p) < d{y,A) and d{p,A) > 0. Let 
X G [p, y). Let a G (0, 00) he arbitrary. Then d{x,p) < d{x, A) — r for any r > 
satisfying 



The proof of Lemma 5.5 is based on the strong triangle inequahty of Clarkson 
[13, Theorem 3]. It is interesting to note that although this inequality was formu- 
lated in the very famous paper [13] of Clarkson, it seems that it has been almost 
totally forgotten, and in fact, despite a comprehensive search we have made in the 
literature, we have found evidences to its existence only in [13] and later in [45]. 



In this section we mention several counterexamples which show that the as- 
sumptions in Theorem 5.1 are essential. 

If the space is not uniformly convex, then the Voronoi cells may not be stable 
as shown in Figures 5-8. Here the setting is point sites in a square in with the 
max norm. 

The positive common lower bound expressed in (1) is necessary even in a square 
in the Euclidean plane. Consider X = [—10,10]^, Pi^/j = {(0,/3)} and P2,/3 = 
{(0,-/3)}, where (3 G [0,1]. As long as /3 > 0, the cell dom(Pi_/3, Pa,/?) is the 
upper half of X. However, if /3 = 0, then dom(Pi 0, P2,o) is X. A more interesting 
example occurs when considering in X the rectangle Pi^/3 = [—a, a] x [—10,-/3] 
and the line segment P2 = [—10,10] x {0}, where a,/3 G [0,1]. If /3 = 0, then 
d{Pifi,P2) = 0, and the cell dom(Pi^o, -P2) contains the rectangle [—a, a] x [0,10]. 
However, if /3 > 0, then this cell does not contain this rectangle. 

The assumption expressed in (2) is essential even in the Euclidean plane with 
two points. Indeed, given /3 > 0, let Pi^p = {{(3,1)} and P2 = {(0,-1)}. Then 
dom(Pi 0, P2) is the upper half space. However, if /3 > 0, then the half space 
dom(Pi^^,P2) contains points {x,y) with y —00. Thus the Hausdorff distance 
between the original cell dom(Pi 0, P2) and the perturbed one dom(Pi^^, P2) is 00, 
so there can be no stability. 



We conclude the paper with the following remarks. First, despite the counterex- 
amples mentioned above, some of the assumptions can be weakened, with some 
caution. For instance, under a certain compactness assumption and a certain geo- 
metric condition that the sites should satisfy, the main result can be generalized 
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Figure 5. Four sites in a square 
in (M2,£oo). The cell of Pi = 
{(0,0)} is displayed. The other 
sites are P2 = {(2,0)},P3 = 
{(-2.0)}. P, = {(0.-2)}. 



Figure 6. Now either P4 is the 
square [-f3, f3] x [-2 - /3, -2 + f3] 
or Pi = {(/?,/?)}, P > arbitrary 
small. The two lower rays have 
disappeared. No stability. 




Figure 7. The full diagram of 
Figure 5. Note the large inter- 
section between cells 1,2, and 3. 
For emphasizing this intersection, 
each cell is represented as a union 
of rays (see Theorem 5.4 for more 
information) and some rays were 
emphasized. 



Figure 8. The full diagram of 
Figure 6 when Pi = {(/3,/3)}. 
Cells 2,3 have been (significantly) 
changed too. 
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to normed spaces which are not uniformly convex. A particular case of these 
assumptions is when the space is m- dimensional with the | ■ |oo norm, the sites 
are positively separated, and no two points of different sites are on a hyperplane 
perpendicular to one of the standard axes. 

Second, an interesting (but not immediate) consequence of the main result is 
that it implies the stability of the (multi-dimensional) volume, namely a small 
change in the sites yields a small change in the volumes of the corresponding 
Voronoi cells. 

Third, it can be shown that the function T defined in Theorem 5.4 also has a 
certain continuity property if the space is uniformly convex and this expresses a 
certain continuity property of the boundary of the cells. 

Fourth, the estimate for A from Remark 5.3 is definitely not optimal and it can 
be improved, but, as simple examples show, A cannot be too much larger and its 
estimate should be taken into account when performing a relevant analysis. There 
is nothing strange in this and the situation is analogous to the familiar case of real 
valued functions. For instance, consider / : M — >■ M defined by /(x) = 0, x < 0, 
/(x) = (l//3)x, X G [0,/3], f{x) = 1, X > f3, where /3 > is given. Although / is 
continuous, a "large" change near x = (of more than /3) will cause a large change 
to /. 
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8. Appendix 1: proof of the main theorem 

This appendix is devoted for the proofs. In the sequel, unless otherwise stated, 
{X, I ■ I) is a normed space; X is a nonempty closed and convex subset of X; 
P, P', A, and A' are nonempty subsets of X; {Pk)k&K and {Pl^k&K are two tuples of 
nonempty subsets of X representing the sites and the perturbed ones respectively; 
it is assumed that the distance between any point x E X and each of the subsets 
P, P', Pfc, P^ (for any k G K) is attained; for each k & K, we let A^ = [Jj^k ^"^^ 
^'k — Ujyfc -Pj j u^it vectors will usually be denoted by 9 or 0. 

The following simple lemma is needed for proving the representation theorem 
below it. 

Lemma 8.1. Let 7^ A C X. Suppose that y,p e X satisfy d{y,p) < d{y,A). 
Then d{x,p) < d{x, A) for any x G [p, y]. 

Proof. Let a G A. Since d{y,p) < d{y,a) and x G [p,y], we have 

d{y,x) + d{x,p) = d{y,p) < d{y,a) < d{y,x) + d{x,a). 
Thus d{x,p) < d{x, a) for each a G A, so d{x,p) < d{x, A). □ 

Theorem 8.2. (The representation theorem) Suppose that for all x & X the 

distance between x and P is attained. Then dom(P, A) is a union of line segments 
starting at the points of P. More precisely, given p & P and a unit vector 9, let 

(3) T{9,p) = sup{t e [0, 00) : p + t9 e X and d{p + t9,p) < d{p + t9,A)}. 
Then 

(4) dom(P,A) = U U [P,P + nO,p)0]. 

peP \e\=i 

WhenT{9,p) = 00, the notation [p,p+T [9, p)9] means the ray {p+t9 : t E [0, 00)}. 

Proof. Given p E P and \9\ = 1, let T{9,p) be defined as in (3). Obviously the set 
which defines T is nonempty and T{9, p) > 0. If T{9, p) = 0, then [p, p+T{9, p)9] C 
dom(P, A). Otherwise, the segment [p,p + s9] is contained in dom(P, A) for any 
s G [0,T{9,p)) by the definition of T{9,p) and Lemma 8.1. If T{9,p) = 00, then 
this simply means that the ray {p + t9 : t E [0, 00)} is contained in dom(P, A), 
so [p,p + T{9,p)9] C dom(P, A) by our notation. Otherwise, T{9,p) is finite. 
To see that p + T{9,p)9 G dom(P, A), we simply use the fact that the function 
t I— !■ d{j) + t9,p) — d{j> + t9,A) is non-positive for all t G [0,T(^,p)), and that 
it is continuous. All the points are in X because it is closed and convex. Thus 
Up6pU|.|=i[P,P + r(^,p)^] C dom(P,A). 

Now let y G dom(P, A). By assumption, d{y,P) = d{y,p) for some p E P. If 
y = p, then y is in the union IJ^gp U|6»|=i fe' 1~'<~T{9, q)9] . Otherwise, let t = d{y, p) 
and 9q = {y — p)/t. Then ^0 ^ : |^| = !}• In order to show that y is in the 
union [jq^p[j^0^^i[q,q + T{9,q)9], it suffices to show that y E [p,p + T{9q,p)9o]. 
Since d{y,p) < d{y,A) and y = p + 19 E X , this is a simple consequence of the 
definition of T{9q,p). □ 
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The following elementary lemma will be useful later. 

Lemma 8.3. Let X be a normed space, 
(a) Given x,y, z E X , x y, z y, 
z — y X — y 
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(5) 



y\ 



\x - y\ 



< min < 2 



\z — x\ 



\x - y\ 



\z — x\ 



(b) Let X C X be nonempty and convex. Let p and y,p^y be two points in X . 
Suppose that the open ball Bh{p,R) = B{p,R) HH, relative to the affine hull 
H spanned by X, is contained in X. If x G {p,y), then Bnix^r) C X for any 
r < mm{d{x,y),d{x,p),Rd{x,y)/{2d{p,y))}. 

Proof, (a) From the triangle inequality, 



z-y 


x-y 




\x - y\{z -y)-\z- y\{x - y) 


\z - y\ 


\x-y\ 




\z — y\\x — y\ 



{\x-y 




z-y\){z 






z-y\{{z-y) - {x-y)) 


2 

< 


\x — z 


\z-y 


1 2 


\x — z\ 




z-y\ 


\x-y 






z-y\ 


x-y\ 




x-y\ 



x\ \z 



y\ is proved 



The inequality \{z — y)/\z — y\ — [x — y)/\x — y\ \ < 2\z 
in a similar way. 

(b) Let z G Bh{x, r). li z = x, then obviously z E X. Now assume that z ^ x. By 
assumption d{y,x) > r and z G B{x,r). Thus y z. Let 9i = {x — y)/\x — y\ = 
{p ~ y)/\p ~ y\ 62 = {z — y)/\z — y\. The point y + d{p, y)92 is in H since 
[y, z] C H. This point is also in B{p, R) since 

\{y + d{p,y)92)-p\ = \{y + d{p,y)92) - {y + d{p,y)9i)\ = d{p,y)\9i - 9^1 

< 2\p — y\\z — x\/\x — y\ < 2\p — y\r /\x — y\ < R 

by part (a) and the choice of r. Therefore y + d{p, y)92 G Bh{p, R) C X. In 
addition, 

d{z, y) < d{z, x) + d{x, y) < r + d{x, y) < d{p, x) + d{x, y) = d{p, y), 
and hence 2 G [y,y + d{j>, y)92\ C X. 



□ 



The following technical proposition is central. 



Proposition 8.4. Suppose that d(P, A) > 0. Let e > be such that e < d{P, A)/ 6 
and suppose that the following conditions hold: 

(6) There exists A G (0, e) such that for each p E P,y E dom(P, A), x E [p, y], 

if d{x,y) = e/2, then d{x,p) < d{x,A) — A. 

(7) ^ 

There exists A' G (0, e) such that for each p' G P', y' G dom(P', A'), x' G [p', y'] 

if d{x',y') = e/2, then d{x',p') < d{x',A') - A'. 
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Suppose that there are two positive numbers M, M' such that 

(8) snp{T{e,p):peP,\e\ = l}<M, snp{T'{e\p) : p' e P',\e'\ = 1} < M', 

where T'(6',p') is defined as in (3), but with A' instead of A. Let €1,63 be positive 
numbers satisfying 

(9) ei + 63(3 + 4M/e) < A/2, ei + 63(3 + AM' /e) < X' /2. 
If 

(10) D{A,A')<ei, D{P,P')<e3, 
then D{dom{P,A),dom{P',A')) < e. 

Proof. We first note that 5e < d{P', A'). Indeed, let p' G P' and a' G A' be given. 
By (10) and (9) there are p & P and a G A such that d{p,p') < 63 < e/2 and 
d{a', a) < ei < e/2, so 

(11) 6e < (i(p, a) < d{p,p') + (i(p', a') + d{a' , a) < e + (i(p', a'), 

which implies that 5e < d{P',A'). We will show that there exists e' G (0, e) such 
that D{dom{P, A),dom{P' , A')) <7 < e. In fact, we will show that we can take 
e' := e — £3. In order to prove that D{dom{P, A),dom{P' , A')) < ?, it suffices 
to show that for any y G dom(P, A) there exists y' G dom(P', A') such that 
d{y, y') < and that for any z' G dom(P', A') there exists z G dom(P, A) such 
that d{z', z) < 7. We will show the first inequality, since the proof of the second 
one is similar using (7), (8), and (9) with A' and M'. 

Let y G dom(P, A). By Theorem 8.2 we can write y = p + 19 for some p G P, 
a unit vector 6, and t G [0,T{e,p)]. By (8) we have t < M. By (10) there is 
p' G P' such that d{p,p') < 63. Suppose first that t < e. Let y' = y. Then y' G X, 
d{y,y') < e — €3, and d{y',p') < d{y',p) + d{p,p') < 2e. In addition, 

(12) 5e < d(P', A') < d{p', A') < d{p', y') + d{y\ A') < 2e + d{y\ A'). 

Hence d{y\p') < 2e < rf(?/', A') and ?/' G dom(P', A'). 

Now assume that e < t. Then t = d{p,y) < d{p,p') + d{p',y) < + \y — p'\. 
Hence < t-es < since ^-63 > t-e/2 > e/2 > 0. Let 6' = (y -p')/\y -p'\ 

and denote €2 = ^e^/e. Then, using Lemma 8.3(a), 



(13) \d-d'\ 



y-p y-p 



\y-p\ \y-p' 



2\p-p'\ 2e3 
\y-p\ e 



Let X = p + {t — 0.5e)6'. Then x G [p, y] and d{x,y) = e/2, so by (6) we have 
d{x,p) < d{x, A) — A. Let y' = p' + {t — 0.5e)6''. Because t — e/2 < \y — p'\ (see the 
discussion above (13)) we have y' G [p\y] C X. In addition, 

%,?/') < d{y,x)+d{x,y') < e/2+(t-0.5e)|e-^'|+rf(p,p') < e/2+Me2+e3 < e-e3 
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by (8) and (9). Recalling that d{y',A) < d{y',A') + D{A,A') holds in general, we 
have 

d{y', p') < d{y', x) + d{x, p) + d{p, p') < Mea + eg + d{x, A) - A + eg 

< Me2 + eg + d{x, y') + d{y\ A) + D{A, A) + eg - A 

< 2Me2 + 3eg + ei + d{y' , A') - \ < d{y', A') - A/2, 

where we used (9) and €2 = 2eg/e in the last inequality. Thus d{y',p') < d{y',A'). 
Consequently, in both cases we have d{y, y') < ? := e — eg and y' G dom(P', A'), 
as required. □ 

Proposition 8.5. Suppose that m.i{d{Pk, Pj) : j,k G K, j ^ k} > 0. Let e > 
be such that e < mf{d{Pk, Pj) '■ j,k e K, j 7^ k}/6. Suppose that the following 
conditions hold: 

(14) 

3 A G (0, e) such that for each k E K,p E Pk, y G dom(Pfc, Ak), and a; G [p, y] 

if d{x,y) = e/2, then d{x,p) < d{x,Ak) — A. 

(15) 

3 A' G (0, e) such that for each k e K,p' e Pi,, y' G dom(P^, A'^), and x' G [p' , y'] 

if d{x\y') = e/2, then d{x',p') < d{x',A',,) - A'. 

Let 

Rk = dom(Pfc, Ak), R'k = dom(P^, A'^). 
Suppose that there are M, M' G (0, 00) such that for all k & K , 
(16) 

sup{Tfc(^,p) : p G P, 1^1 = 1} < M, snp{U{9',p') : p' G P', \9'\ = 1} < M' , 

where Tk{0,p) and Tl,{9',p') are defined as in (3), but with Ak and A'^ instead of 
A. Let 64 be a positive number satisfying 

(17) 4(1 + M/e)e4 < A/2, 4(1 + M7e)e4 < A72. 
// 

(18) D{Pk,Pl)<e, ykeK, 
then D{Rk, R'k) < e for each k E K . 

Proof. From the condition e < d{Pk, Pj)/Q for each k ^ j it follows that e < 
^(-Pfc, Ujy/c In addition, the assumption D{Pk,Pk) < for aS\. k E K 
implies D(|J^g^ Pj, |J^g^ P/) < for any I C K, and in particular this is true for 
/ = K\{k}, i.e., D{Ak,A'i,) < e^. Indeed, we only have to observe that given 
y G IJjgj Pi, there exists j G / such that y E Pj, so 

d{y,[jP:) < diy,P;) < DiP„P;) < e4, 

and the inequality follows. Let eg := e4 and ei := e4 + r, where r is defined as 
the minimum of the two positive numbers (A/2 — 4e4(l + M/e))/4 and (A'/2 — 
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4e4(l + M7e))/4. Then ei + 63(8 + 4M/e) < A/2 and ei + €3(8 + 4M7e) < AV2. 
In addition D(Pfc,P^) < €3 and -D(y4fc,v4^) < £4 < ei for each k E K. Thus all 
the conditions of Proposition 8.4 are satisfied for each k E K with P = and 
A = Ak. Hence D{Rk, Ri) < e. □ 

The following lemma provides a simple sufficient condition for (8) to hold. 

Lemma 8.6. Suppose that the following condition holds: 

(19) 3p G (0, 00) such that Vx G X the open ball B{x,p) intersects A. 
Then 

(a) snp{T{e,p) ■.\e\ = l,peP}< p. 

(b) d{x, A) < p and d{x, A') < p + D{A, A') for each x e X. 

(c) sup{T'{e',p') : l^'l = l,p' e P'} <p + D{A,A'). 

In the same way, if A^ = Uj^^ Pj for each k E K and the following condition 
holds: 

(20) 3p G (0, 00) such that for all /c G -ftT and for all x G X 

the open ball B{x,p) intersects Ak, 

then all the above claims remain true with Tk,Pk,Ak instead ofT,P,A. 

Proof, (a) Let p G P be arbitrary and let 6 be an arbitrary unit vector. Let t = p 
and X = p + t6. If x ^ X then T{6,p) < t, since otherwise there is some s > t 
such that p + s9 E X and d{p + s9, p) < d{p + s6, A) by the definition of T (see 
(3)). Hence x G [p,p + s6] C X, a contradiction. Now assume that x E X. By 
(19) there exists a E A such that d{x,a) < p = d{x,p), so d{x,A) < d{x,p). 
By the definition of T and Lemma 8.1 we conclude that T{6,p) < p in this 
case too. 

(b) Let X E X. By (19) we know that d{x, A) < p. Since d{x, A') < d{x, A) + 
D{A, A') holds in general, we have d{x, A') < p + D{A, A'). 

(c) If D{A, A') = 00, then the assertion holds trivially. Assume now that D{A, A') < 
00. By part (b) the open ball B{x, p + D{A, A')) intersects A' for each x E X 
and the proof continues as in part (a) with p + D(A, A') instead of p. 

The proofs in the case where (20) holds are the same as above. □ 

The next lemma is the key step in establishing (6), (7), (14), and (15) in the 
case of uniformly convex normed spaces. It improves upon Lemma 8.1 under the 
additional assumptions that {X, | ■ |) is uniformly convex and d{p, A) > 0. Its proof 
is based on the following definition and on a special case of the forgotten strong 
triangle inequality of Clarkson [18, Theorem 8]. See page 6 for the definition of 
uniformly convex spaces and the meaning of the function 6. Clarkson's theorem, 
which is formulated in [18] for uniformly convex Banach spaces, has nothing to do 
with completeness and it remains true without this assumption. 

Definition 8.7. Given two non-vanishing vectors x,y E X, the angle (or Clark- 
son's angle, or the normed angle) a{x,y) between them is the distance between 
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their directions, i.e., it is defined by 

fI \y\ 

Theorem 8.8. (Clarkson) Let xi, X2 be two non-vanishing vectors in a uniformly 
convex normed space (X, | ■ |). //xi + X2 7^ 0, then 

(21) \xi + X2I < + \x2\ - 25{ai)\xi\ - 25{a2)\x2\, 
where ai = a{xi, xi + X2), I = 1,2. 

Lemma 8.9. Let (X, | ■ |) be a uniformly convex normed space and let (!) A C X . 
Suppose that y,p E X satisfy d(y,p) < d{y, A) and d{p, A) > 0. Let x G [p, y). Let 
a G (0, 00) be arbitrary. Then d{x,p) < d{x, A) — r for any r > satisfying 

.oo^ ^ • / A) f d{p,A) 

(22) r < mm < (J, ,d{y,x)d ' 



10 ' \10{d{x,A) + a + d{y,x))^ 

Proof. The assertion is obvious if either p = y ot p = x, so assume p ^ y and 
p ^ X from now on. Let a G A satisfy d{x, a) < d{x, A) + r. We distinguish 
between two cases. Suppose first that the angle a{x — y,a — y) is greater or equal 
to d{p, A)/ {10 {d{x, A) + a + d{y,x))). By the assumption, the strong triangle 
inequality (21) and the fact that 6 is increasing, 

d{y, x) + d{x,p) = d{y,p) < d{y, A) < d{y, a) 

< d{y, x) + d{x, a) — 26{a{x — y,a — y))d{y, x) — 26{a{a — x, a — y))d{x, a) 

< d{y, x) + d{x, A) + r - 2r - 0. 

Hence d{x,p) + r < d{x,A) as required. We note that all the angles are well 
defined, since < d{x,p) < d{x,a) by Lemma 8.1, and < d{y,p) < d{y,a), < 
d{y, x) by assumption. In addition, the expression in (22) is well defined, since 
d{p, A) < d{p, x) + d{x, A) < 2d{x, A) by Lemma 8.1, so the argument inside 6 is 
in the interval [0,0.2] C [0,2]. In addition, the minimum in (22) is positive since 
all the numbers inside 6 are positive and the space is uniformly convex. 

Assume now that a{x — y,a — y) < d{p, A)/{10{d{x, A) + a + d{y, x))). Since x 
is between p and y, the unit vector 9 = {p — y)/\p — y\ satisfies p = y + d{y,p)6 
and X = y + d{y,x)9. Let = (a — y)/\a — y\. Let z = y + d{y,a)9. Since 
d{y,p) < d{y, a), it follows that p is between z and x (and maybe p = z). Because 
a{x—y, a—y) = d{9, 0) and r < a, it follows from the assumption on a{x — y, a — y) 
that 

A ^fa ^^d{p,A){d{y,x) + d{x,a)) d{p,A) 

diz.a) = d(9,(p)d(y,a) < — -— < . 

^ ' ' V ,v^; vy> ; 10{d{x,A) + a + d{y,x)) 10 

Thus 

d{p, A) < dip, a) < dip, z) + diz, a) < dip, z) + dip, A)/10, 
so 9(i(p, A)/ 10 < dip, z). Since p is between z and x, the choice of r in (22) implies 

dix. A) + Adip, A)/10 > dix. A) + r > dix, a) > rf(x, z) - dia, z) 

= dix, p) + dip, z) - dia, z) > dix, p) + Mip, A) /lO - dip. A) /lO. 
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Therefore d{x,p) + r < d{x, A) in this case too. □ 

We are now ready for formulating and proving the main stabihty result (denoted 
by Theorem 5.1 in the main body of the text). 

Theorem 8.10. Let (X, | ■ |) be a uniformly convex normed space, and let X 'O X 
be closed and convex. Let {Pk)k€K, {Pk)k€K be two given tuples of subsets of X 
with the property that the distance between each x E X and each Pk, P^ is attained. 
Suppose that 

f] ■= inf{d(Pfc, P,) : J, /c e K,3 ^ fc} > 0, 
and let A^ = [Jj^j^Pj, A'j^ = [Jj^f,Pj for each k E K. Suppose also that (20) 
holds. For each k E K let Rk = dom{Pk, Ak), R'f, = dom(P^, A'^) be, respectively, 
the Voronoi cells associated with the original site Pk and the perturbed one P^. 
Then for each e G (0,?7/6) there there exists A > 0, namely A = min{Ce^, 0.5e} 
for any positive C satisfying 



16(p + 5?7/12) Vl2p + 5r/^ 

such that if the inequality D{Pk, P(,) < A holds for each k E K, then D{Rk, R'f.) < e 
for each k E K . 

Proof. Let e E (0,?7/6) be given. We will show that all the conditions needed in 
Proposition 8.5 are satisfied. Let k E K he given. As in the first lines of the proof 
of Proposition 8.5, we have D{Ak,A'f}j < sup{D{Pj, Pj) : j ^ k} < A, and hence 
D{Ak,A'f^) < e < 7]/Q. Let M = p and M' = p + r]/Q where p is from (20). By 
Lemma 8.6(a) and Lemma 8.6(c) (with P = Pk, P' = P^ A = Ak, A' = A'J we 
obtain (16). 

Let p' E PI, be given. We claim that 5?7/6 < d{p',A'iJ. Indeed, let a' E A'j^ be 
arbitrary. Since D{Pk,Pk) < A < e/2 and D{Ak,A',^) < e/2 + /3 for any (3 > 0, 
there are p E Pk and a E Ak such that d{p,p') < e/2 and d{a',a) < e/2 + /3. 
Therefore 

V < d{p, Ak) < d{p, a) < d{p,p') + d{p', a') + d{a' , a) < e + (3 + d{p', a'), 

Since a' and /3 were arbitrary and since brj/G < rj — e we have 5r]/Q < d{p',A'j^) 
as claimed. Since 6e < 77 we also conclude that 5e < d{p',A'i^). In addition, 
0.5e < 2e < 4d{p,Ak)/10 and 0.5e < 2e < 4d(p',A'J/10 for any p E Pk and 
P' e Pi 

Let 0" = 77/6. Then obviously 0.5e < e < a. By Lemma 8.6(b) we know that 
d{x,Ak) < p and d{x,A'i) < p + r]/6 for each x E X. As a result, from the 
monotonicity of 5, < 5 < 1, 5?7/6 < mm{d{p,Ak),d{p',A'i,)}, and from Lemma 
8.9 it follows that 

(23) X' = X = 0.5e-6 (—^ 

Vl2p + 5r/, 

satisfy (14) and (15) (the minimum in the expression in the right hand side of 
(22) is the third element and both A and A' are smaller than this element). Let 
e4 = A. Then < Ce^ < eA'/(8(M' + r//4)) < A7(8(l + M'/e)) and hence (17) 
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is satisfied. As a result, we conclude from Proposition 8.5 that D{Rk,R'i^) < e, as 
required. □ 

From Theorem 8.10 we see that A = O(e^). However, under the (practical) 
additional assumption that the sites are strictly contained in the interior of the 
world X, we can obtain a better bound, namely A = 0(e). The proof is based on 
the following assertions. We recall that dX is the boundary of X relative to the 
affine hull spanned by X and that the distance from a point to the empty set is 
infinity. 

Proposition 8.11. Suppose that d{P,A) > and let e > satisfy e < d{P,A)/6. 
Suppose that d{P,dX) > and d{P',dX) > 0. Assume that (6) and (7) hold and 
that there are two positive numbers M, M' such that (8) hold. Suppose that ei and 
€3 are two positive numbers such that the following inequalities hold: 

ei + Seg < min{A/2,AV2}, 

(24) . [ d{P, dX)e d{P',dX)e \ 

^ \ 4 max{M, M'} ' 4 max{M, M'}}' 

If 

(25) D{A,A')<e,, D{P,P')<e,, 
then D{dom{P,A),dom{P',A')) < e. 

Proof. The proof is in the same spirit as the proof of Proposition 8.4 and many 
details are similar and hence they have been omitted here. First, exactly as in 
the first lines of proof there we have 5e < d{P', A'). As explained there, it suffices 
to show that for any y G dom(P, A) there exists y' G dom(P', A') such that 
diy,y') < e := e - 63. 

Let y G dom(P, A). By Theorem 8.2 we can write y = p + t9 for some p G P, a 
unit vector 6, and t G [0, T{e,p)]. By (8) we have t<M. By (25) there is p' G P' 
such that d{p,p') < €3. Suppose first that t < e. Let y' = y. Then y' G dom(P', A') 
exactly as in the proof of Proposition 8.4, and obviously d{y, y') < e — e^. 

Now assume that e <t. Let x = p+{t—0.5e)9. Then x G [p, y] and d{x, y) = e/2, 
so by (6) we have d{x,p) < d{x, A) - X. Let 9' = 9 and y' = p' + {t - e/2)9'. Then 
d{x,y') = d{p,p') < €3, and hence y' ^ X hj Lemma 8.3(b) and the second 
inequality in (24) (in this connection note that since A < e it follows that €3 < 
0.5e < {d{x,y),d{x,p)}). In addition, 

d{y, y') < d{y, x) + d{x, y') = e/2 + d{p, p') < e/2 + 63 < e - €3 

by (24). Recalhng that d{y' , A) < d{y', A') + D{A, A') holds in general, we have 

diy', p') < d{y', x) + d{x, p) + d{p, p') < 63 + rf(x, A) - A + 63 
< 2e3 + d{x, y') + d{y' , A') + D{A, A') - A 

< 3e3 + ei + d{y\ A') - A < d{y' , A) - A/2, 

where we used the first inequality in (24) in the last inequality above. Thus 
d{y'.iP') < d{y',A') and then y' G dom(P',yl'). Therefore, in both cases we have 
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d{y, y') < € = €-€3 and y' G dom(P', A'), and hence D{dom{P, A), dom(P', A')) < 
e. □ 

Proposition 8.12. Suppose that inf{d{Pk, Pj) '■ j,k & K, j ^ k} > and let 

e > be such that e < mi{d{Pk, Pj) '■ j,k ^ K,k ^ j}/6- Suppose also that 
d{\J^^j^Pk,dX) > and d{\Jf^^j^ Pl^,dX) > 0. Suppose also that (14) and (15) 
hold. Let 

Rk = dom(Pfc, Ak), R'k = dom(P^, A^). 
Suppose that there are M, M' G (0, 00) such that for all k & K , (16) hold. Let €4 
he a positive number satisfying 

, ^ . [A d{[j,,KPk^dX)e ^(U,^ P^, gX)6 ' 
^ ^ ' \8'8' 4max{M,M'} ' 4 max{ M,M'} 

// 

(27) D{Pu,Pl)<e, WkeK, 

then D{Rk, R'f,) < e for each k & K. 

Proof. From the condition e < d{Pk, Pj)/Q for each k ^ j it follows that e < 
'^(-Pfc) Ujyfc In addition, as in the proof of Proposition 8.5 the assumption 
P)(Pfe,Pfc) < €4 for a\\ k e K implies that D{Ak,A'jJ < 64 for all k e K. Let 
€3 := €4 and ei := + r, where r is defined to be the minimum of the two positive 
numbers (A/2 - 4e4)/4 and {\'/2 - 4e4)/4. Then (24) is satisfied for P = Pfc for 
each k & K because of (26). In addition D{Pk,P[.) < 63 and D{Ak,A'i^) < £4 < ei 
for each k G K. Thus all the conditions of Proposition 8.11 are satisfied for each 
keK. Hence D{Rk, P'J < e. □ 

Theorem 8.13. Let {X, | ■ |) be a uniformly convex normed space, and let X (1 X 
be closed and convex. Let {Pk)k&K, {Pk)keK be two given tuples of subsets of X with 
the property that the distance between each x & X and each Pk and Pj. is attained. 
Suppose that d{\J^^j^ Pk.dX) > and d{\J^^j^ P^^dX) > 0. For each keK 
let Ak = UjVfc^i, ^fc = Uj^kPj let Rk = dom{Pk, Ak), R'k = dom(P^, A^) be, 
respectively, the Voronoi cells associated with the original site Pk and the perturbed 
one PI.. Suppose that 

7] := inf{rf(P,, P,) -.j^ke K,j ^k}>0, 

and that (20) holds. Then for each e > satisfying e < min{?7/6, 8-(i(|J^,g^ Pk, dX)} 
there exists A > 0, namely A = Ce for any C > satisfying 



16 Vl2p + 5r/y' 8(p + r//6) 

such that if the inequality D{Pk, P'k) < A holds for each k G K , then D{Rk, R'k) < e 
for each keK. 

Proof. The proof is almost identical to the proof of Theorem 8.10 with the excep- 
tion of a small verification which is explained in the next paragraph. The values 
of A and A' are the same those given in (23) and one verifies that with M = p, 
M' = p + rj/Q, and €4 = A the conditions of Proposition 8.12 are satisfied. 
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The small verification that mentioned above and should be done is that the 
inequality < ci(Ufcg^ P^, (9X)e/(4max{M, M'}) holds. This is true for the fol- 
lowing reason: if dX = 0, then the right hand side in this inequality is infinity. 
Otherwise, let b G dX and p' G [j^eK ^'k given. Then p' G for some k ^ K. 
Since 5 < 1 it follows that C < 1/16 and hence < e/16 < 0.5d(\J^^j^ Pk,dX). 
Because D^P^, P^) < it follows that d{p',p) < for some p G Pk and therefore 

d{ [j Pfc, dX) < d{p, h) < d{p,p') + d{p', 6) < €4 + d{p', b) 

<0.5rf(|J Pk,dX) + d{p',b) 

k£K 

Since p' and b were arbitrary it follows that 0.5(i(|J^gj^ P^, dX) < rf(IJfce_K- dX). 
Recalling that 4max{M, M'} = 4(p + 77/6), we conclude that 

^ ^ 0.5-c^(U,^Pfe,gX)6 ^ d{\J,^^Pl,dX)e 
4(p + r//6) - 4(p + r//6) 

□ 

9. Appendix 2: A continuity property of T 

This appendix is devoted to a discussion on the continuity of T(-, p) in uniformly 
convex spaces which was mentioned briefly in Section 7. This result is related to 
the main result and the approach used for deriving it is similar. 

The following definition plays an important role. 

Definition 9.1. Let X he a closed and convex subset of a normed space. Let 
p e X. Let Qp := {6 : \6\ = 1, p + t6 e X for some t > 0} be the set of all 
directions such that rays emanating from p in these directions intersect X not only 
in p. Let 9 G 6p. Let L{9) G (0, 00] he the length of the line segment generated 
from the intersection of X and the ray emanating from p in the direction of 9. 
The point p is said to have the emanation property ( or to satisfy the emanation 
condition) in the direction of 9 if for each e > there exists /3 > such that for 
any cf) G Op, if \(j) — 9\ < P, then the intersection of X and the ray emanating from 
p in the direction of (j) is a line segment of length at least L{9) —e. In other words, 
L{(j)) > L{9) — e. The point p is said to have the emanation property if it has the 
emanation property in the direction of every 9 G Op. A subset C of X is said to 
have the emanation property if each p ^ C has the emanation property. 

As an illustration, in the following examples each point p in the set X satisfies 
the emanation condition. First example: X is any bounded closed convex set 
and p G X is an arbitrary point in the interior of X relative to the affine hull 
spanned by X; second example: the boundary of the bounded closed and convex 
X is strictly convex (if a 7^ 6 are two points in the boundary, then the open line 
segment (a, b) is contained in the interior of X relative to the affine hull spanned 
by X) and p G X is arbitrary. Any ball in a strictly convex space has a strictly 
convex boundary; third example: A is a cube (of any finite dimension) and p G A 
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is arbitrary; fourth example: X is a closed linear subspace and p G X is arbitrary. 
The proof of the third example is by direct computation. The proof of the first and 
the second examples is based on Lemma 9.2 below and, e.g., on Lemma 8.3(b). 

Before formulating Lemma 9.2 we note that the reason we consider only Qp 
and not the whole unit sphere is that in general it may happen that L{(j)) = 
for unit vectors (p arbitrary close to 6 (consider e.g., X = [0,1]^ in the Eu- 
clidean plane, p = (0,0), 6 = (1,0), and = (^1 — l/ra^, — 1/n), n G N or 
X = [0, 1]2 X {0} in with the Euclidean norm, p = (0,0,0), 6 = (1,0,0), and 
(f)n = (a/1 — 0,—1/n), n G N) and hence the condition L{(j)) > L{6) — e does 
not hold. Intuitively, rays in the directions of (p ^ Qp go out X immediately and 
they do not affect the "events" inside X. 

Lemma 9.2. Let X be any closed and convex subset of a normed space, let p G X , 
and let 6 & Qp be given. Assume that L{6) < oo and that any point in the segment 
{p,p + L{9)9) is in the interior of X relative to the affine hull H spanned by X. 
Then p satisfies the emanation condition in the direction of 9. 

Proof. L{9) > since 9 G Qp. Let e G (0, L{9)) be given. Let z = p + {L{9) - e)9. 
Then z is well defined and it is in the interior of the segment [p,p + L{9)9] and 
hence there exists p > such that the intersection of the ball B{z, p) and the affine 
hull H is contained in X. Let r = p/{L{9) — e). For each G 6p, if — 9\ < r, 
then |(p+(L(^)-e)0)-^| = (L(^) -e)|0-^| < p. Thus the point p+ (L(^) -e)0, 
which is in H since (p G Qp, is also in the ball B{z,p) and hence in X. Thus 
L((/)) > L{9) - e by the definition of L(0). □ 

Example 9.3. The emanation condition does not hold in general. Consider the 
Hilbert space £2- Let (e„,)^]^ be the standard basis. Let yi = Ci and for each n > 1 
let yn = ei/2 + e„/n. Let A = {— Ci}. Let X be the closed convex hull generated 
by A [j{yn : n = 1, 2, . . .}. Let p = and let 9n = yn/\yn\ for each n. Then p does 
not satisfy the emanation condition in the direction of 9i. Indeed, lim„_^oo^n = ^1 
but L{9n) = v/0.25 + 1/^2 < 0.99 = L{9i) - 0.01 for each n>l. The subset X is 
in fact compact since the sequence {yn)'^=i converges. 

Using the definition of the emanation condition, we can prove the continuity of 
T{-,p). The proof is based on several lemmas. 

Lemma 9.4. Let X be a convex subset of a normed space. Let p & X and 7^ 
A G X be given. Suppose that d{p,A) > 0. Let e G {0,d{p, A)/6) and assume that 
the following conditions hold: 

(28) There exists A G (0, e) such that for each y G dom(p, A), x G [p,y], 

if d{x, y) = e/2, then d{x,p) < d{x, A) — A, 

(29) there exists M G (0, 00) such that snp{T{9,p) : \9\ = 1} < M. 

Let P be a positive number satisfying (3 < X/ {AM). Let yi = p + ti9i, i = 1,2 where 
ti>t2 + 0.5e, t2 > 0, and 9i e Qp, i = 1,2. If yi G dom(p, A) and \9i - 6*21 < /3, 
then d{y2,p) < d{y2,A). 
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Proof. Suppose first that ti < e. Then 

6e < d{p, A) < d{p, y2) + d{y2, A) < e + d{y2, A) 



since t2 <ti < e. Therefore d{y2,p) < e < 5e < d{y2, A) as claimed. 

Now consider the case where e < ti. Let x = p + {ti — e/2)6i and z = p + 
{ti — e/2)^2- Note that z and also y2 are not necessarily in X, but this is not 
important. By (29) we have ti < M. Hence d{x,z) = (ti - e/2)\9i - 6*21 < M/3. 
Since d{x,yi) = e/2, yi G dom{p,A), and x G [p,yi], it follows from (28) and the 
choice of (3 that 



Since t2 < ti — 0.5e it follows that y2 G [p, z]. As a result, Lemma 8.1 implies that 



Lemma 9.5. Let X be a closed and convex subset of a normed space. Let A ^ X . 
Let p E X and let 9 E Qp. Suppose that there exists M G (0, oo) such that 
T{(f),p) G [0,M] for all (f) G 6p, where T{(j),p) is defined with respect to A in (3). 
Then for any e G (0, oo) there exists (3 > such that for all cj) G 0p, if \0 — < (3, 
then T{(j),p) < T{e,p) + e. 

Proof. Suppose that this is not true. Then for some e G (0, oo) there exists a 
sequence (0„)^]^ of elements in 9p such that 9 = lim^^oo 0n but T{(j)n,p) > 
T{9,p) + e. By assumption T(0„,p) G [T{9,p) + e, M + e] for all n. Therefore there 
exists a subsequence {(j)nj)fLi such that t := \imj^ooT{(j)nj,p) exists and satisfies 
te [T{9,p)+e,M+e]. Hence T(0„^,p)0„^ t9 and Urn j^oc{p+T{(f)nj,p)(pn,) = P+ 
t9. Because X is closed it follows that p+t9 G X. By Theorem 8.2 (with P = {p}) 
we know that p + T(0„^.,p)0„^. G dom.{p,A) for each j. Therefore the inequality 
d{p + T{(f)n^)(j)nj,p) - dip + T{(j)npP)<PnpA) < holds for each j. Thus, from the 
continuity of the distance function it follows that d{p + t9,p) — d{p + t9, A) < 0. 
As a result p + t9 E dom(p. A). But T{9,p) < t since T{9,p) <t — e, and this is a 
contradiction to the definition of T{9,p) (see (3)). □ 

Theorem 9.6. Let X be a closed and convex subset of a uniformly convex normed 
space and let A be a subset of X . Let p E X and suppose that t] := d{p,A) > 0. 
Suppose that (19) holds and that p has the emanation property. Then for each 
e G (0, 77/6) and each 9 E Qp there exists A E (0, 00) such that for each G Qp, if 
\9-4>\< A, then \T{9,p) - T{(P,p)\ < e. 

Proof. Let cr = 77/6 and 



Then A satisfies (28). Indeed, let x and y be such that y E dom{p,A), x E [p,y], 
and d{x,y) = 0.5e. We have d{x,A) < p by Lemma 8.6(b) since (19) holds. In 
addition, < A < 0.5e < a < Arj/lO and A is not greater than the minimum in 
(22). From Lemma 8.9 it follows that d{x,p) < d{x,A) — A. 



d{z, p) < d{z, x) + d{x, p) < M(3 + d{x, A) - \ 



<Mf3 + d{x, z) + d{z, A) - A < d{z, A) - A/2. 



d{y2,p) < d{y2,A). 



□ 
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From Lemma 9.5 there exists /3i > such that for all cf) G 0p, if \9 — (f)\ < (3i, 
then T{(j),p) < T{9,p) +e. Since p has the emanation property there exists /92 > 
such that for all G Op, if \6 - < /32, then L{(j)) > L{9) - e. Let 

A = min{/3i,/32,A/(4p)}. 

We claim that for each G 6p, if l^' — 0| < A, then \T{6,p) — T(0,p)| < e. Indeed, 
let G 6p satisfy |^ — 0| < A. Then \9 — (f)\ < so Lemma 9.5 implies that 
T(0,p)<T(^,p) + e. 

Assume for contradiction that T{9,p) > T{(f),p) + e. Let ti = T{9,p), yi = 
p + ti9, t2 = ti — e, and 1/2 = P + ^20- In particular ^2 > T{(f),p). Because 
1^-01 < it follows that L(0) > L(^) - e > r(^,j9) - e = ^2- Therefore the 
point 1/2, which is on the ray emanating from p in the direction of 0, is in X. 
By Lemma 8.6(a) (with P = {p}) we know that (29) is satisfied with M = p. 
In addition, since \9 — (f)\ < A/(4p) and since yi E X and d{yi,p) < d{yi,A) 
(by Theorem 8.2), it follows from Lemma 9.4 that d{y2,p) < d{y2,A). By the 
definition of T{(j),p) we therefore have ^2 < T{(f),p), a contradiction. □ 

Remark 9.7. The conditions mentioned in Theorem 9.6 are necessary for the 
continuity of T mentioned there. Indeed, the uniform convexity of the norm is 
needed as can be seen in the case described in Figure 5 where p = (0, 0) and 
A = {(0, —2), (2, 0), (—2, 0)}. The discontinuity of T(-,p) is at the four unit vectors 
9 = (±1,±1). The emanation property is necessary since in Example 9.3 it can 
be seen that T{-,p), p = 0, is discontinuous at 6* = ^1. As for the assumption 
d{p,A) > 0, consider X = [—1,1]^ in the Euclidean plane where p = (0,0) and 
A = [—1, 1] X {0}. Then T{9,p) is discontinuous at 6^ = (0, 1) because T{9,p) = 1 
but T{(f),p) = for any of the form = {1/n, ^1 — l/n^), n G N. As for the 
boundedness condition (19), let X = with the Euclidean norm, p = (0,-1), 
and A = {(0, 1)}. Then T{-,p) is not continuous at 6* = (1, 0), because T{9,p) = 00 
but T{(j),p) < 00 for any = (a/1 — l/n^, l/n), n G N. 



